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1 Introduction 



In the past year, there has been substantial progress in the study of scattering amphtudes 
in three-dimensional Chern-Simons-matter theories [1-5]. Much of this work has focused on 
the = 6 superconformal Chern-Simons theory developed by Aharony, Bergman, Jafferis 
and Maldacena (ABJM), since in a certain limit, this theory is dual to type IIA string theory 
on AdS4 X CP^, and therefore provides a new example of the AdS/CFT correspondence 
[6]. The AdS4/CFT3 correspondence has many features in common with the AdS5/CFT4 
correspondence, which relates J\f — ^ super Yang-Mills theory (sYM) to type IIB string 
theory on AdSs x [7]. This suggests that many of the symmetries, structures, and 
dualities exhibited by the planar amplitudes of = 4 sYM may also present in the ABJM 
theory. 

The first hint of this possibility was the discovery that the four- and six-point tree-level 
amplitudes of the ABJM theory have Yangian symmetry [2]. In [8], it was demonstrated 
that the Yangian symmetry of = 4 sYM is equivalent to superconformal plus dual 
superconformal symmetry. Dual superconformal symmetry is hidden from the point of view 
of the action, and can only be seen by studying on-shell amplitudes [9, 10]. The equivalence 
of Yangian symmetry with superconformal plus dual superconformal symmetry when acting 
on amplitudes was also demonstrated for the ABJM theory in [5]. Furthermore, an integral 
formula similar to the Grassmannian integral formula of A/" = 4 sYM [11] was proposed for 
the ABJM theory in [4]. This formula has Yangian symmetry and is conjectured to generate 
all tree-level amplitudes. This was verified for the four-point amplitude and claimed to 
be verified for the six-point amplitude [4]. In order to demonstrate dual superconformal 
symmetry and test the Grassmannian integral formula for amplitudes with more than six 
external legs, we need an efficient way to compute higher-point on-shell amplitudes in the 
ABJM theory. Unfortunately, explicit Feynman diagram calculations become tedious for 
amplitudes with more than six external legs. 

The computation of amplitudes can be drastically simplified using linear identities, such 
as supersymmetric Ward identities, which relate amplitudes with the same number of legs 
and helicity structure in four dimensions [12, 13], or off-shell nonlinear recursion relations 
such as the Berends and Giele recursion relations [14], which relate higher-point amplitudes 
to products of lower-point amplitudes. In recent years, twistor-inspired methods [15] have 
led to new representations of tree-level amplitudes such as the MHV vertex expansion [16] 
and the Britto, Feng, Cachazo and Witten (BCFW) [17, 18] recursion relations, both of 
which give a systematic procedure for constructing higher-point amplitudes from lower- 
point amplitudes. 

The BCFW recursion relation is derived by analytically continuing the momenta into 
the complex plane. In particular, one shifts the momenta of two external particles by a 
complex parameter. The analytic continuation is such that the amplitudes become rational 
functions in the complex plane, whose only singularities arise from propagators becoming 
on-shell. The recursion relation then follows from the fact that a rational function can be 
reconstructed from the residues of its poles, provided that the function vanishes at infinity. 
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The residues in this case are simply products of lower-point on-shell amplitudes. 

One of the crucial ingredients of the BCFW formalism is the complex deformation, 
which is defined such that (super-) momentum conservation and on-shell properties of the 
external fields are preserved. As we describe in section 3, however, the standard deforma- 
tion used in the BCFW approach is only valid for D > 4. In this paper, we will define 
a complex deformation such that all properties of the on-shell amplitudes are preserved 
in three dimensions. To construct a recursion relation from this deformation, one must 
demonstrate that the amplitudes vanish when the parameter of the deformation goes to 
infinity, which is not always true even for D > 4. We argue that the superamplitudes 
of the ABJM theory indeed have vanishing large-z behavior. We confirm the validity of 
our recursion relation by using it to reproduce six-point amplitudes computed in ref. [2] 
using Feynman diagrams and an eight-point amplitude computed from the Grassmannian 
integral formula given in ref. [4]. 

While an on-shell recursion relation provides an efficient way to construct higher-point 
amplitudes from lower-point amplitudes, it also serves as a convenient tool for analyzing 
the symmetries of amplitudes. If one can demonstrate that the recursion relation pre- 
serves symmetries of the lowest-point tree-level amplitude, it follows by induction that the 
symmetry holds for all tree-level amplitudes. Indeed, the proof of dual superconformal sym- 
metry of four- [19], six- [20], and ten-dimensional [21] maximal sYM planar amplitudes was 
based on the fact that the BCFW recursion formula preserves dual superconformal symme- 
try. Equipped with a similar recursion relation in three dimensions, and the fact that the 
four- and six-point ABJM tree amplitudes are dual superconformal invariant [1, 5], we will 
prove that the recursion relation also preserves the dual symmetry, and hence all tree-level 
amplitudes of the ABJM theory have dual superconformal symmetry. Using generalized 
unitarity methods [22-25], we further extend this symmetry to the cut constructible part 
of the loop amplitude. 

Dual superconformal symmetry in = 4 sYM is related to the fact that type IIB 
string theory on AdSs x is self-dual after T-dualizing certain directions corresponding 
to bosonic and fermionic isometrics [26-28]. The dual superconformal symmetry of all 
tree-level amplitudes in the ABJM theory therefore strongly suggests that type IIA string 
theory on AdS4 xCP^ is self-dual after T-dualizing certain bosonic and fermionic isometrics, 
although various attempts to demonstrate this have encountered singularities [29-33]. 

For = 4 sYM, the combination of ordinary and dual superconformal symmetry is 
sufficient to fix the tree-level amplitudes [34]. Since the Grassmannian formula for the 
ABJM theory has both of these symmetries, the fact that the ABJM tree-level amplitudes 
have dual conformal symmetry suggests that they are indeed generated by the Grassman- 
nian formula. 

The structure of this paper is as follows. In section 2, we review the spinor-helicity 
formalism for three-dimensions. In section 3, we explain the difficulties of extending the 
BCFW recursion relation to three-dimensions and derive a new recursion relation for three- 
dimensional theories. The novel feature of this recursion relation is that it requires deform- 
ing the momentum spinors of two external particles in a nonlinear way. In section 4, we use 
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background field methods to argue that the superamphtudes of the ABJM theory vanish 
when the deformation parameter is taken to infinity. In section 5, we review the Grass- 
mannian integral formula for the ABJM theory and show that it reproduces the six-point 
superamplitude computed in ref. [2] using Feynman diagrams. We also use it to compute 
an eight-point amplitude for the first time. In section 6, we use our recursion relation to 
reproduce the six-point and eight-point amplitudes computed in section 5. In section 7, 
we use our recursion relation to prove for the first time that all tree-level amplitudes of 
the ABJM theory have dual superconformal symmetry. We extend this result to loop-level 
using generalized unitarity methods in section 8. We conclude with some discussion in 
section 9. In appendices A and B, we review our conventions and derive the Feynman rules 
for the ABJM theory in the background field approach. In appendix C, we describe an al- 
ternative recursion relation for three-dimensions which requires deforming the momentum 
spinors of four external particles in a linear way. We also present a superconformally co- 
variant form of the recursion relation derived in section 3. Appendix D provide additional 
computational details for various results in the paper. 

Note added: Immediately after the first version of this paper was submitted to the 
e-print archive, we were informed of the paper [54] which also studies the scattering am- 
plitudes of the ABJM theory and its mass deformation. 

2 Three dimensional spinor helicity and on-shell superspace 

Here we give a brief introduction to the three dimensional spinor helicity formalism [1-3]. 
A vector in three dimensions, when written in bi-spinor form, is given by a symmetric 2x2 
matrix. Each spinor transforms under the three-dimensional Lorentz group Spin(l,2) = 
SL(2,R). A null momentum in three dimensions can be written in bi-spinor form using a 
single spinor: 

= A^A^. (2.1) 

Our conventions for spinors and gamma matrices are summarized in appendix A. Here, 
we simply note that for real p^, the gamma matrices are such that p^^ is also real. The 
spinor A*^ is real for outgoing (j)o < 0) particles and purely imaginary for incoming (j)o > 0) 
particles. 

The spinors are contracted via the SL(2, R) metric 

(ii)^A«A,„ = Afe„^Af, (2.2) 

where ei2 = — e^^ = 1. The vector and spinor Lorentz invariants are related by 

2pi-Pi = -(uf . (2.3) 

In later sections, we will make use of the "A-parity" symmetry A —A. It does not 
change the momentum p^^ ^ but flips the sign of the fermion wave- function. As a result, it 
acts on amplitudes by 

^(Ai, . . . , -Ai, . . . , A„) = (-l)^^^(Ai, . . . , Ai, . . . , A„) , (2.4) 
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where Fi is the fermion number of the particle on leg i. We note that for Chern-Simons 
matter theory, only amplitudes with even number of external states are non-vanishing. 
This is because for odd number of states, there will always be at least one Chern-Simons 
gauge field on the external line. Since the gauge field does not carry any physical degree 
of freedom, the amplitude vanishes. 

The on-shell superspace for the AB JM theory is built upon three fermionic coordinates 
77^ in addition to A*^ [2], which transform as a 3 under the U(3) subgroup of the S0(6) 
R-symmetry group. The particle/anti-particle superfields take the form 

$ = 04 + ^^^^ + leijKV^ij^cj)^ + l^uWv^V^^A , 

$ = ^4 ^ + leuWr^'i^"" + leijWr^'r^''4>A • (2.5) 

The SO (6) R-symmetry acting on the J\f — Q supercharges are realized by 77^ and their 
conjugates = d/drj^ through the Clifford algebra, 

{ri'Xj}-S'j (/,J=1,2,3). (2.6) 

It is customary [2, 4] to introduce a collective notation A = (A; 77). The color ordered 
amplitude, following the convention of ref. [4], is then a function of A^ and we choose to 
associate ^odd/ewen to multiplet. Then, the supersymmetric extension of eq.(2.4) is 

^(Ai, • • • , -Ai, • • • , A„) = (-l)M(Ai, . . . , Ai, • • • , A„) . (2.7) 

The generators of the superconformal symmetry come in three types: 

The invariance of the amplitudes under the first type of generators will be manifest from 
the notation. For the second type, we will use the notation 

_ ^al ^ ^a^I^ ^IJ ^ ^ ^2.9) 

The invariance under p^^ and q"^^ will be imposed by the super-momentum conserving 
delta functions 

6Hp)S\Q) with P^J2pf^ Q^J2lf- (2.10) 

i i 

A similar delta function does not exist for r^^. In ref. [2], it was shown that the r^^ 
invariance introduces a coset 0(2A: — 4)/U(A: — 2) for the 2A:-point amplitude. As we will 
review in section 5, the same coset appears in the Grassmannian integral formula of ref. [4]. 

3 Recursion relation 
3.1 Momentum shift 

In the L)-dimensional BCFW formalism (with > 3), one introduces a complex parameter 
z into the amplitude by shifting the momenta of two of the external legs [18, 35, 36]. For 
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legs i and j, the shift is given by 

Pi^Pi + zq, pj^pj-zq, (3.1) 

for some vector q. For the external legs to remain on shell, the shift vector q is required to 
satisfy 

q • Pi = q • Pj = = ' (3.2) 

Although these constraints admit nontrivial solutions in D > 3 dimensions, for D — the 
only solution is g = 0. To see this, recall that a null vector in three dimensions can be 
written as a bi-spinor product A'^A^. Since the spinor is two dimensional, without loss of 
generality one can write the spinor for the vector q as 

\q = a\i + h\j . (3.3) 

One then immediately sees that the requirements in eq.(3.2) imply a = 6 = 0. 

The analysis above can be translated into the spinor language. Under a general two-line 
shift A A(z), momentum conservation implies that 

A^A^ + A^-A^- = Xi{z)Xi{z) + Xj{z)Xj{z) . (3.4) 

Assuming that the shift is a linear transformation on (A^, Aj), one can write 

Then momentum conservation requires that 

R^{z)R{z) = /, (3.6) 

i.e. the shift is an element of S0(2,C). ^ 

Note that the BCFW deformation in four dimensions can also be formulated in terms 
of a shift matrix. In four dimensions, the momentum conservation can be written as 

n n 

5^p, = 5^A,A, = 0, (3.7) 

i=i i=i 

which is invariant under an arbitrary GL(n, C) transformation, 

Xi^MiXj, Xi^XjiM-^i. (3.8) 

The BCFW shift (Aj — )■ Aj + zXj, Xj — ^ Aj — zXi) can then be regarded as a particular 
element of GL(2, C) C GL(n, C): 



Xi{z)\ ^ 1 z\ Xi 
X,{z) [Ol [X, 



(3.9) 



^ Since we are dealing with a continuous deformation, we cannot use 0(2, C) including orientation reversal. 
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Starting from the three dimensional momentum conservation and following the same logic, 
one naturally arrives at the S0(2,C) shift as above. 

Let us now return to the case of three-dimensions. If we assume that R{z) in eq.(3.5) 
is linear in then it can be parameterized as 

\ a2iz 1 + a22Z J 

It is not difficult, however, to see that eq.(3.6) constrains all of a's to be zero. A similar 
analysis shows that eq.(3.6) also constrains the three-line shift to be trivial, i.e. if we take 
R{z) to be an S0(3,C) matrix that depends linearly on z eq.(3.6), constrains it to be the 
identity. The first non-trivial solution which depends linearly on z appears at four-lines. 
Multi-line shifts, when compared to two-line shifts, have the disadvantage of introducing a 
larger number of diagrams for each recursion. Furthermore, for the purpose of analyzing 
the perseverance of symmetries throughout the recursion, one would prefer as few shifted 
variables as possible. However, since multi-line shifts generally have better large z behavior, 
which may be necessary for lower supersymmetric theories, we give a detail account of the 
four-line shift in appendix C. 

Alternatively, we could introduce a mass deformation. In that case, the spinor helicity 
formalism is essentially four-dimensional, so one can define a linear two-line shift. For the 
purpose of proving dual superconformal symmetry of the amplitudes, however, it would be 
better to avoid introducing a mass-deformation. 

In order to construct a two-line shift without introducing a mass deformation, we must 
relax the assumption that R{z) is linear in. z. In the next section, we will construct such 
an S0(2, C) matrix and use it to derive a recursion relation for tree-level amplitudes. 

3.2 Derivation of the recursion relation 

We parameterize the R{z) matrix by 

/ z-\-z~-^ z—z~^ \ 



m = \j.-^ .j-^: • (3.11) 



2i 



It acts on two reference external legs, say, the i'th and the j'th. Using the cyclic symmetry, 
we can set z = 1 and the j = / with 2 < / < A: + 1 without loss of generality. Since we 
are interested in a recursion relation for superamplitudes, we also need to consider super- 
momentum conservation, 

n n 

Y.^f-Y.^Ui-^- (3-12) 

The conservation of both momentum and super-momentum can be maintained if we also 
deform the 77's [19]. Thus we define the super-shift as: 



- m ? , '^^ri - Riz) "M , (3.13) 
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or collectively we have 

(3.14) 




Note that, since R^R = 1, the super-shift preserves not only the (super-)momentum but 
also all the other superconformal generators, 

EA.^. EA.A.. Es^- (3.15) 

i i i 

We will denote the tree-level superamplitude for given external particles and their 
momenta as 

^(Ai, • • • , A2k) = A{Ai, ■■■ , A2fc)^'(P). (3.16) 

After the deformation, we can regard it as a function of z. Following the original proof 
of the BCFW recursion relation [17], we express the undeformed amplitude at z = 1 as a 
contour integral along a small circle around 2: = 1, 

A{z^l)^^<[ (3.17) 

The on-shell tree- level superamplitude has the following factorization property [2], 
^|,.^o - /^3 /.(A/,---,A.,A,)A^(iA;,A.^,,---,A.-,) ^ ^^^^^^^^ ^ 

^ J Pf 



(3.18) 



in the limit 0, where pf = pj + • • • + pj is the momentum of the factorization 

channel /. The left/right sub-amplitudes Al,Ar are on-shell at = 0. The variable 
Af = {XJ;r]^) is allocated to the internal propagator, where Xf is determined by the 
momentum conservation. The relation (3.18) holds up to an unspecified overall constant. 

Now we apply the 50(2) deformation R{z) to eq.(3.18) and plug the resulting A{z) 
into the integrand of eq.(3.17). Then, we change the contour of eq.(3.17) to encircle the 
poles arising from the deformed internal propagator l/pf{z)'^. As we will show in the next 
section, the integrand vanishes sufficiently fast at z ^ oc and z ^ 0.^ Thus we obtain 

^(. = l) = --^> W A>M ^ +6 \ J^ AUz-,rj)A^{z;iv) _ ^3_^g^ 




pj{z) 

This is the three-dimensional J\f = 6 supersymmetric recursion relation. The sum over all 
possible channels is denoted by X^j, and (ztz* jr) are zeroes ofpj{z). The recursion relation 
for component amplitudes can be obtained by applying the corresponding 77-integration to 
both sides of eq.(3.19). We will see explicit examples in the next section. 



^ A{l/z) is the same as A{z) up to a sign because exchanging z with l/z in eq.(3.11) acts as an orientation 
reversal. Sufficiently fast vanishing means A{z 00) ^ 0{l/z) which implies A(z ^ 0) ^ ^(^)- 
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In our deformation, the on-shell condition j3j(z) = takes the fohowing form 

pj(z) = afz-^ + ^/ + CfZ^ = 0, (3.20) 

with aj, 6/, Cf independent of z. To see this, let the deformed momenta be pi and pi. They 
are in the fohowing form, 

pi = AiAi = ^{pi +pi) + z^q + z~^q, 

pi = = ^(pi ^pi) - z^q - z~^q, (3.21) 

where q and q are given by 

= i(Ai + iAO"(Ai + iAO^ = ^(Ai - iAO"(Ai - iA;)''. (3.22) 

Using the cychcity, let us choose / and J to satisfy 1 < / < / < J, then Pf{z) is given by 
Pf{z) — Pi ^ • • • + Pi{z) + • • • pj. From this, we obtain 

af = -2g • {pf - pi) , hf = (p/ + pi) • {pf - pi) Cf = -2q • (p/ - pi) , (3.23) 

after some manipulation using q^q— \{pi—pi)^ Q' {Pi^Pl) — Q' {Pi ~^Pl) — 0^ etc. Using 
eq.(3.23), one can explicitly write down the zeros {izj^ ^,±^2/} eq.(3.20): 



' (-2,/) } - I 4,.(p,-pO / ' ^'-''^ 

where we used a variation of Schouten's identity, 

{r\p\s)'^ = (r|p|r)(5|p|5) +p^{rs)'^^ (for each r, 5). (3.25) 

Let us rewrite the sum over residues of eq.(3.19) so that the factorization limit is 
transparent: 

A{z^l)^J2 [ d'7i\{H{zlj,z;j)AL{zlj;v)An{zlj;i7i) + {z^j o zlf)) , (3.26) 

^ Pf 



f 



where the function H{a, b) is defined by 



The second equality defines the function /i(a, b). We used the fact that 

Al{-z)Ar{-z) = {-iy+^ALiz)ARiz) , 

which follows from eq.(2.7). Since h(a^b)AL(a)Aji(a) is invariant under both a —a and 
b —6, it is a function of o? ^b^ . The numerator in eq.(3.26), 

^K/,^2*j)^L«/)^i?«/) - h{zlj,zlj)AL{zlj)AR{zlj), 

is antisymmetric in ^ ^ z^ j. Given that it is a function of (z^ j)^ and (2^2 j)^, it should be 
divisible by (z^ j)^ — (2^2 j)^- Moreover, since the quotient must now be a symmetric function 
under z\ ^ ^ 2^2/5 it can be written as a rational function containing (^1 j)^ + (^2/)^ ~ 
—bf/cf and (zj^ j)^(2:2 j)^ = ^f/^f- Thus one concludes that the final result is free from 
any square-root factors that {±^^^,±2:2/} contain. 
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4 Asymptotic behavior of the superamphtude 

In this section, we study the asymptotic behavior of the deformed superamphtude under 
the deformation in eq.(3.13). The change of integration contour from eq.(3.17) to eq.(3.19) 
is justified only when A{z) vanishes as z ^ oc or 0, so that there is no pole at z = oo and 
z = 0. We will focus on the behavior as z ^ oc, as the behavior for z ^ can be derived 
in a similar fashion; see footnote 2. 

Under a BCFW deformation, the large-z behavior of a gluon amplitude in = 4 sYM 
depends on the helicity of the shifted legs. Although certain helicity configurations lead 
to bad large-z behavior, when all the component amplitudes are packaged into a super- 
amplitude, the entire superamphtude vanishes at large z if the supershift is appropriately 
defined. 

Here, we work in the opposite direction. In particular, we first identify component 
amplitudes that share the same large-z behavior as the superamphtude. Then we show that 
these component amplitudes behave as ^ as z ^ oo, which implies that the superamphtude 
also vanishes at large z. 

Choosing a component amplitude We need to clarify what it means to say that a 
component amplitude has the same large-z behavior as the superamphtude. In principle, 
the large-z behavior of the superamphtude receives contributions from different component 
amplitudes, each of which has different large-z behavior, so it is not obvious that one can 
make a direct connection between the large-z behavior of particular component amplitudes 
and the superamphtude. 

To answer this question, consider shifting legs i and j of the superamphtude, and 
expand the superamphtude in (r]i^r]j): 

A = + Af''\{ + Af'^^vj + ■■■ + A^'''\vi)Hvj? . (4.1) 

Each sub-amplitude ^(^'^) depends on all A's and all 77's except {rji^rjj). After the super- 
shift, (A^, Aj) {Ai{z), Aj{z)), the superamphtude becomes 

A{z) = + A\^'^\z)fji {z) + Af'^\z)fijiz) + ■■■ 

= ^(O'O) (z) + 4^'°) iz)v! + 4°'') iz)vj + ---. (4.2) 

On the first line, the z-dependence of A^''^{z) is entirely due to the shift (Ai,Aj) 
{Xi{z), Xj{z)). Thus, these sub-amplitudes contain the large z behavior of the compo- 
nent amplitudes under the bosonic shift. On the second line, the expansion variables are 
the undeformed (77^, 77^). The relation between A^^'^\z) and A^^'^\z) follows from the 
shift property of {r]i,r]j). For example, and ^(^'^^ are singlets in the sense that 

^(0,0) ^ ^(0,0) ^(3,3) ^ ^(3,3)^ ^j^g Q^j^gj. hand, we have 

P'°^^^)Vm.)-K'°^^^)1 (4 3) 

[a^'^'\z))-^^'^ 1^^(0-1) (z) J • 
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Figure 1. Diagrammatic representation of a background field computation. 



Now, the key point is that all A^^'^\z) share the same large-z behavior even though 
A^^'^\z) do not. To see this, recall that the super-shift preserves the super-momentum 
conservation condition, or equivalently the supersymmetry generator Q^^ = qf^^) . On 
the other hand, the individual supersymmetry generators qf^^ can relate different sub- 
amplitudes A^'^''^\z), which implies the supersymmetric Ward identity. In particular, 
applying super-momentum conservation and collecting terms linear in 77/ and independent 
of rjj, we find 

J2qf^A{z)^0 ^ = (no sum in/), (4.4) 

where = "^i^ij if^- Since the prefactor on the right-hand side has no 2-dependence, 
^(O'O) and A^^'^^ share the same large-z behavior. Applying the same method iteratively, 
we can show that all A^'^''^\z) share the same large-z behavior, which then defines the 
large-z behavior of the full superamplitude A(z). 

Note that ^(^'^^ only includes component amplitudes with legs i,j corresponding to 
either ^4{Xi) or '04(Aj), since these are the lowest weighted components in the superfield 
expansion in eq.(2.5). Since A^^'^\z) = A^^'^\z), their large-z behavior is identical. Hence, 
the large-z behavior of the superamplitude is the same as the large-z behavior of component 
amplitudes of the form (• • • (^^4(1) or '04(i)^ • • • or '04(j)^ • • • ). 

4.1 Background field method and naive large-z behavior 

We will now demonstrate that component amplitudes of the form (• • • '04(i) • • • 04(j) • • • ) 
vanish as 1/z in the large- 2: limit. To simplify the analysis, we will take the shifted 
legs to be adjacent. This fixes the amplitudes under consideration to have the form 

('^4(1) 04(0)- Note that the shifted legs do not have to be adjacent in order to 

have good large-z behavior. In the next section, we will use the Grassmannian integral 
formula to show that non- adjacent shifts result in even better large-z behavior. Another 
consequence of the adjacent shift is that A{—z) = —A{z) as follows from the "A-parity" 
(2.7). We will see that several apparently 0{z^) terms all cancel out. 

In the large-z limit, it is convenient to analyze these amplitudes using a background 
field formulation which describes hard particles scattering through a soft background [35]. 
Since we are shifting two external legs, the amplitudes under consideration reduce to dia- 
grams of the form depicted in Fig.l, where the crosses represent background fields. 
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We proceed by separating each field into a background contribution plus a quantum 
fluctuation: 

^ + a^, ^ + 0^ ¥ ^¥ + t/j^. 

There are similar equations for the hermitian conjugates of the matter fields.^ While the 
background fields are solutions to the classical equations of motion, we take the fluctuations 
to describe fields whose momentum has been shifted and therefore approaches z^g as z ^ 
oc, where q is the term in the deformation of external momenta pi and pi, 

q-^ = 1 (A? + zXf) (A^ + zAf ) . (4.5) 

We define A^ = ^(A^ + zAf ). 

Since all the vertices in Fig.l involve at most two fluctuation fields, it is sufficient to 
expand the action to quadratic order in fluctuations for the purpose of analyzing diagrams 
of this type. Schematically, such terms take the form 

a'^A, a2$2^ a(f>^A, A^cj?, Axp'^, xpa"^, ^^^^ <^^*^, 

a'^(a^</.)$, A^'{d^<P)<P, a''{d^^)<p. (4.6) 

Note that terms linear in the fluctuations vanish by the equations of motion for the back- 
ground fields. Let us analyze their leading large-z behavior individually. The first row 
contains the kinetic terms for the quantum fields. The corresponding propagators can be 
shown to behave as (1,^,1) respectively. Hence, only the scalar propagators improve the 
large-z behavior. The second and third rows contain interaction terms. Note that terms 
in the second row do not contain any derivatives, so their vertices do not contribute any 
powers of z. 

We now focus on the third row. The a^{d^^)(j) vertex does not introduce z-dependence 
since the derivative is on a background field. On the other hand, the a^(5^0)$ vertex 
behaves as . Since a gauge field fluctuation cannot appear as an external leg in the 
diagrams that we consider, any vertex involving gauge field fluctuations must be connected 
to a matter vertex via a vector propagator. The combination of the a^(5^(/))$ vertex and 
a single vector propagator behaves as: 

z^oo [p + z'^qY 2z'^q-p 

Hence, this combination approaches a constant as z ^ oc. Finally, the large-z dependence 
of the A^{dij^(j))(j) vertex can be set to zero if gauge-fix the background so that 

g-A = 0. (4.8) 



^Note that the ABJM theory has two gauge fields A, A for the U(iV) x U(A^) gauge group. In this 
section, we will not distinguish between them, since this does not play an important role in the analysis. 
For simplicity, we will also neglect color-indices and R-symmetry indices when they are not essential to the 
discussion. 
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There are some circumstances for which this gauge choice cannot be implemented. As 
explained in ref. [35], this problem will arise for diagrams in which the shifted external 
lines interact with a single background gauge field via a three-point contact term. This 
situation does not arise in our analysis, however, because we take the shifted external lines 
to be a scalar and a fermion, which cannot interact with a single background gauge field 
via a three-point contact term. 

From the above analysis, we see that the propagators and vertices in eq.(4.6) will at 
worst give a constant in the large-z limit. Combined with the fact that there is an external 
wavefunction for the fermion in ('04(— z^g) • • • ^4(2:^^)), which behaves as lim^^oo '04(^^^) — 
zXq, a naive analysis would imply that the amplitudes can have at most an 0{z) divergence 
at large z. On the other hand, since the scalar propagator introduces a factor of ^, 
diagrams containing at least one scalar propagator are 0{l/z) in the large-z limit, so 
the only diagrams that we need to worry about are those which do not contain scalar 
propagators. As we will see in the next subsection, however, the large- 2: behavior of these 
diagrams is 0{l/z) as well. 

4.2 Improved large-z behavior 

To see the effect of the external wave function, we analyze all possible ways of connecting 
the external ip4{z) line to the rest of the diagram. We focus on diagrams without scalar 
propagators, and we will show that each case behaves as ^ asymptotically. Since from the 
previous discussion we've concluded that the building blocks of the amplitude behave at 
worst as a constant, the ^ behavior of the sub diagram then implies the entire diagram can 
behave at worst as ^. 

Fermion propagator: If the external hard fermion is connected to the rest of the di- 
agram via there will be a fermion propagator. The asymptotic of the two 
scenarios is 



z^oo {p + z'^qy 2z'^q-p z' 

z^oo + z^g)^ zz'^q • p z 

where we used the notation {ip{p)\ and to represent the external wave function and 

the background field respectively. For both cases we have used the fact that the leading 
large z contribution in the external wave function is ip4{pi) zXq^ and q-A{k) = q^^\(3 — 
0. 

Gluon propagator: If the external hard fermion is connected to the rest of the diagram 
via il^gl^ ^ there will be a gluon propagator. Since the gluon will be required to turn into 
a scalar at some point along the line, the other side of the propagator will be required to 
connect to a^(5^0)$, a^^^A^, e^^^a^a^^A^ or a^l>$a^. 



-13- 




X 



(a) 



(b) 



Figure 2. Gluon propagator connecting to scalar vertex cancels against contact vertex in the 
leading order. 




Figure 3. Gluon propagator connecting to fermion vertex. 



1. We first consider the case where the gluon propagator directly connects to the 
a^(5^0)$ vertex. Such diagrams generically take the form of Fig. 2 (a). ^ We note 
that for every Fig. 2 (a) there will be complimentary diagram using the four-point 
contact vertex ^^'00$. as in Fig. 2(b). The asymptotic behavior is of order z for both 
diagrams, and they cancel each other. The details of this cancelation is given in 
appendix D. However, it is not difficult to be convinced of this result once one notes 
that if the background fields are replaced with on-shell external lines, it is simply the 
four-point amplitude which we know vanishes as ^ at large z. 

2. The gluon propagator can also connect to a fermion vertex, and change to scalar later 
down the line. This correspond to the diagram in Fig. 3. The second fermion-vertex 
will have to connect to a fermion propagator. Using the Feynman rules given in 
appendix B, the two vertex two propagator combination gives 

(.10) 

The leading large z behavior takes the form: 

z{q\a^^{hf-^^{k2)a^ (4.11) 
q-ki q- {ki + ks) 

which vanishes using the identity eq.(A.7). 

3. Next we consider the gluon propagator connecting to e^^pa^ay . This corresponds 
to the diagram in Fig. 4 and gives 



^Note that we have denoted the scalar to be 04, even though this is not always the case. The justification 
is that if the scalar is of a different R-index, then one needs to connect to other scalar vertices to change 
the R-index. This will introduce scalar propagators. 
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Figure 4. Gluon propagator connecting to three gluon vertex. 




\ 




(X) 

Figure 5. Gluon propagator connecting to scalar propagator via quartic vertex. 



-g' {m)\<^'m f':'}^'^':i^ (e'^nA.ih f ^^^^ . (4.12) 



(Pl + h) {k3+P2 

For the leading z contribution, pi = P2 = z^q, and the computation reduces to 



-zg^ {q\ ^f^^^ki f^-^'^^''^'"M^^^'--^l' + o{l/z) 

2q ■ kiq ■ fcs 

= ^^'^ff^f^ ('^l <^^^{k,)e^sq' + 0{\lz). (4.13) 

where we've used — Since the order z piece is proportional to 

q • A{k2)'> it vanishes for our background field gauge, and hence the diagrams start at 
1 

z ' 

4. For the case where the gluon propagator is connected to the vertex a^^^A^^ as 
indicated in Fig. 5, there are two types of vertices depending on the type of gauge 
fields on the four-point vertex. The computation for the two are the same, so we 
present the result for the first diagram: 

(4.14) 

In the large-z limit this becomes 

-zg^ (g| a^^k. f^^^f^y^^ + OiXjz) 



2q-ki 



^q- A{k2){q\^{ki) 



^-zg^^-^f^f^ + Oil/z). (4.15) 

where we have used eq.(A.7). Thus again the leading z pieces vanish due to the 
background field gauge. 
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Figure 6. Gluon propagator connecting to another gluon propagator via quartic vertex. 



5. Finally, we consider the gluon propagator connected to the vertex a^$$a^. This is 
shown in Fig. 6 and gives 

The leading large z behavior takes the form: 

^ /"^y ^ ^ . (4.17) 
Contracting the Levi-Cevita tensors, this term vanishes since = j^Xq = 0. 

Thus we see that for all cases, due to the fact that leading large z behavior of the 
external fermion wave function is proportional to Ag, the asymptotic behavior of amplitudes 
with adjacent (^4,^4 as shifted lines are improved to K Since these amplitudes share the 
same large z behavior as the superamplitude, our recursion relation is applicable. In the 
next section, the same large z behavior will be reproduced by shifting the Grassmannian 
integral formula. Finally, note that for the limit z ^ 0, the same analysis would give that 
the amplitude behaves as z and hence vanishes as well. 

5 Grassmannian integral 

In this section, we digress to discuss the matrix integral formula for the tree level amplitudes 
of planar AB JM theory proposed in ref . [4] : 

^ ^~ J vol[GL(A:)] MiM2---Mk_iMk * 

The integration variable C is a (A: x 2k) matrix. The dot products denote (C • C^)mn = 
CmiCnii {C • A)^ = Cmi^i- The z-th minor Mi of C is defined by 

Mi = e'^^'"'^^C^i(i)C'^2(^+i) • • • Cmk{i^k-i) • (5.2) 

In ref. [4] , this formula was shown to satisfy the same cyclic symmetry and superconformal 
symmetry as the tree-level (2/c)-point superamplitude. It was also shown to reproduce the 
known 4-point superamplitude and satisfy Yangian invariance. 

After a brief review of its defining properties, we will use the formula to revisit the 
asymptotic z behavior of the previous section and then to compute some amplitudes up to 
8-point. For both purposes, the formula reveals remarkably simple structures that are not 
visible in a Feynman diagram analysis. 



■'2k 
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Orthogonal Grassmannian The integral (5.1) should be understood as a contour in- 
tegral on the moduli space of rank fc, {k x 2k) matrices Cmi (m = 1, . . . , A;; i = 1, . . . , 2k) 
subject to the constraint {C • C^)mn = CmiCni = and the equivalence relation C ^ gC 
with g G GL(A;). This moduli space is known as the orthogonal Grassmannian 0G(A:,2A;). 
The constraint and GL(A:) gauge symmetry determines the dimension of OG(A:, 2k) as 

dime [OG(fc, 2k)] = 2fc2 - fc2 - M^±i^ = kik-Xl 

Here, we give an informal introduction to the orthogonal Grassmannian.^ For more formal 
treatments, see, for instance, refs. [37, 38]. 

In the mathematics literature, OG(/c, 2k) is defined to be the moduli space of maximal 
isotropic subspace of C^^. We assume that C^^ is equipped with a non-degenerate quadratic 
form (y, W) for V^W ^ C^^ and choose to work in a basis of C^^ such that {V, W) = 
ESi ^i^i' A subspace X C C^^ is called isotropic if {V, W) ^ for ah V,W ^ X. The 
maximal dimension of such X is known to be k. If we interpret the {k x 2k) matrix C 
as a collection of k vectors in C^^, the constraint C • — and the equivalence relation 
C ^ hC, h ^ GL(A:) dictate that the k vectors form a basis of a maximally isotropic 
subspace X. 

Alternatively, OG(A:, 2k) can be defined to be the coset space: ^ 

0G{k,2k) = 0{2k)/V{k). (5.4) 

Geometrically, it is the moduli space of complex structures of M?^ compatible with the 
standard metric. To see this, we introduce complex orthonormal frames in M?^ satisfying 

ds^ = dx'dx' = {Ermdx'){E^dx^) ^ ® , E^ = {Erm)"" , 

^mii^j) ^ij 5 EjniE^ 25^ , EmiEjii = = E^ E^ . (5.5) 

In terms of the orthonormal frames, the complex structure J^j can be written as 

J^j = S^'Jij , i dx' A dx^ ^E"^ AEm- (5.6) 

The complex structure is compatible with the metric by construction. It is also clear that 
the moduli space of J^j is 0(2A:)/U(/c), where 0{2k) and U(/c) acts on Ej^i from the right 
and left, respectively. 

The two definitions of OG(fc, 2k) can be related to each other as follows. To go from 
the coset to the C matrix, we simply take Cmi — ^mi- To define the inverse map, note 
that C • is a positive definite hermitian matrix. By diagonalizing it, we can write 



CC^ = U^DU 




^We thank Hoil Kim for discussions on basic facts about Grassmannian. 

^As we will see shortly, for the integral formula to work, it is important to include both of the two 
disconnected S0(2/c) components of 0(2k). 
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So, E = ^/2/DUC defines the inverse map. It is possible to show that the map between 
(the equivalence classes of) C and E is one-to-one, but we will not give the details here. 

We can discuss coordinate patches of OG(A:, 2k) in the same way as for the ordinary 
Grassmannian. We choose k out of 2k columns of C and fill them with the (k x k) identity 
matrix Ij^. There are (2A:)!/(A:!)^ different patches. For instance, we can have 

C=(^4c), C-C^^O =^ cc^ + 4 = 0. (5.8) 

The (fc X k) matrix c is (zbz times) a complex-valued orthogonal matrix. Thus each coor- 
dinate patch is isomorphic to 0(A:,C). 

After taking account of the bosonic delta function 6{C • X) and the overall momentum 
conservation, the net number of integration variables of eq.(5.1) becomes [4] 

M^-a + 3^'^--^f-^'. (5.9) 

The delta function constraint Cmi^f = means that (A^^^, A^^^) regarded as two vectors 
in C^^ are contained in the maximal isotropic subspace spanned by C. Following steps 
similar to those of the four dimensional case [11], it is possible to do a partial gauge 
fixing by choosing two other vectors in C^^ such that the remaining integration variables 
parameterize OG{k — 2,2k — 4) whose dimension is given by eq.(5.9). ^ 

5.1 Asymptotic z revisited 

If we consider R{z) as an element of 0{2k, C) by the natural embedding 0(2, C) C 0(2A:, C), 
it acts on the k x 2k matrix integration variable C from the right. We can use this fact to 
rewrite the deformed Grassmannian integral as 

(K( f '^"'"'"^ (5^(C-C'^)(5^^l3fe(C-i?(z)-A) 
L2k{A(z)) - J ^^^(.L(fc)] Mi(C) • • • Mk{C) 

^kx2k^ 5^(C'-C'^)52fc|3fc(C'.A) 



~ J Yo\[Gh{k)]Mi{C ■ R-^{z))---Mk{C ■ R-^{z))' ^ ~ ' 

-I 



vol[GL(A:)] Mi{C{z))---Mk{C{z)) 



{C{z) = C-R-\z)) . (5.10) 



On the second line, we used the GL(fc) gauge invariance. On the third line, we removed the 
tilde to avoid clutter. Note that all z-dependence has been absorbed into the minors. This 
is in accord with the discussion in section 4 that the full superamplitude has a well-defined 
large z-behavior. 

As in section 3.1, suppose the legs 1 and / are deformed. From the z-dependence of 
the minors, we now argue that the asymptotic behavior of the superamplitude is 

A{z) 0(z-(^-i)) , (2 < / < A:) . (5.11) 



The importance of this coset in constructing R-symmetry invariants was noted in ref. [2]. 
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To see this, first note that 

C± = Ci±iQ =^ c^{z) = ^C+ + ^C., Ci{z) = ^.C+-^C.. (5.12) 
Here, Ci {1 < i < 2k) denote the z-th column in C. Mi is the determinant of 

(diiz),C2,--- ,Ciiz),--- ,Ck), 

thus Ml is independent of z. Miyi are also independent of z since the determinant includes 
neither of Ci(z). On the other hand, each M2<i<i{z) are linear in z at the leading 

order as z ^ oc, thereby confirming the z~^^~^'^ behavior claimed in eq.(5.11). 

For the case / = /c + 1, there is a subtlety in finding the large z behavior from the 
Grassmannian integral. Recall that the integral (5.10) is actually a contour integral. After 
solving the bosonic delta function constraints, we can write the contour integral schemati- 
cally as 

where r parameterize the {k — 2){k — 3)/2 holomorphic coordinates of OG(A: — 2,2k — 4). 
The function J(t) denotes the Jacobian factor arising from the delta function constraints. 
From eq.(5.12), it is easy to see that (r dependence is suppressed for simplicity) 

Miiz) = |(det(C+,C2,...,Cfc) + 0(z-2)) =zmi + 0(z-i) , 
M2iz) = ^ (det(C2, ...,Ck, C+) + Oiz-^)) = zm2 + 0{z-') , 
Miiz)^^^{det{Ci,...,Ck,C+,...,Ci+k-i) + Oiz-^)) , for 2 < i < fc. (5.14) 

One peculiarity of the I = k + 1 case is that the zeros of the two minors Mi(r; z), M2(r; z) 
tend to collapse at large z. Since Ci{z) and iCi{z) both approaches (2;/2)(C+ + 0{z~'^)), 
if we denote the zeros of Mi(r; z) by r*, then 

ri* = r2* + 0(^-2). (5.15) 

At T^, all minors Mj^i are 0{z) at large z, except M2 which is 0{z~^)\ 

M2(t*; z) = zm2(T2* + 0{z-^)) + O^z'^) = 0{z-^). (5.16) 

Thus residue at r^, denoted by {1}, asymptotically behaves as 

{1} 0{z-^ X ) = 0(z-(^-2)). (5.17) 

z~^ x z . . .z 

First z~^ factor comes from the residue of 1/Mi at r^. The residue at t|, {2}, shares the 
same property. Since at r* {i > 2) all minors Mj^i are 0{z), we conclude that 

{i} = 0{z-^), 2<i<k. (5.18) 
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Thus depending on the choice of contour, the asymptotic behavior of superampHtude could 
be or 0{z-'') for I ^ k + 1. In any case, it vanishes as z goes to infinity. For 

/c = 2, 3, there is no integration, so the superamphtudes behave as 0{z~^) at large z. 

If the Grassmannian integral is indeed capable of producing all tree-level amplitudes, 
the above observation shows that the supersymmetric recursion relation should work for 
any choice of reference legs. This is consistent with (and stronger than) the analysis of 
section 4 based on Feynman diagrams. 

We can also consider a non-supersymmetric recursion relation by deforming A only and 
leaving 77 intact. The asymptotic behavior of the amplitudes can become worse, since the 
fermionic delta function can contribute up to z^ order to the numerator of C2k{z)' Similar 
supersymmetric improvements have been noted when applying the BCFW deformation to 
four dimensional = 4 sYM [19], and AA = 8 supergravity [39]. 

Finally, eq.(5.11) tells us that the asymptotic behavior of the amplitude becomes better 
when we deform two legs further apart. This tendency is also observed in four dimensions, 
which can be used to derive the generalized BCJ relations [40] as demonstrated in [41]. 

5.2 Computing amplitudes up to 8-point 

4-point amplitude The 4-point superamplitude was first computed by Feynman dia- 
grams in ref. [1]. In our convention, the amplitude is given by 

(14)(34) • ^^-^^^ 

It is rather trivial to reproduce this amplitude from the Grassmannian integral. One can 
fix the GL(2) gauge by setting, for instance, 

C = h ' . (5.20) 

\cai C43 ly 

The integral reduces to finding the zeroes of the bosonic delta function 5{C • A). 

Inserting c^^ to the remaining delta functions and the minors, one finds 

Combining eq.(5.21) and eq.(5.22), we recover the superamplitude eq.(5.19). For later 
convenience, we list some explicit examples of component amplitudes here. 

(<^(Ai)</.(A2)<^(A3)0(A4)) = ^^^^ , (^(Ai)</>(A2)^(A3)^(A4)) = > (5-23) 

(<^(Ai)^(A2)^(A3)0(A4)) = -j^^^ , (^(Ai),/.(A2)<^(A3)V'(A4)) = ^^^y . (5.24) 
Throughout this section, {(l),(f>,ip,xj)) without superscripts wih denote (0^, (^4, -04, '^*^). 
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5.2.1 6-point amplitude 

Generality The counting (5.9) for the number of integration variables shows that there 
is no integral to do other than solving the delta function constraints. But, since 0G(1, 2) = 
0{2)/U{l) = Z2, the amplitude is expected to be a sum of two terms. In fact, the general 
structure of the 6-point superamplitude was analyzed in ref. [2]. There, based on the 
invariance under the SO (6) i?-symmetry, it was shown to take the following general form, 

Ae = S^P)5\Q) {5'{a+)U{X) + 5'{a.)f.{X)) . (5.25) 

Here, 6^ {a) = ^ejjKO^^o^'^o^^ and = a\.[X)ril is some linear combination of rfl that is 
linearly independent of the Q-directions. The two terms are related to each other by the 
"A-parity" symmetry. 

As noted in ref. [2], there is some ambiguity in the precise definition of a±. From our 
perspective, the ambiguity is nothing but the GL(/c) gauge symmetry of the Grassman- 
nian integral. It was also shown in ref. [2] that superconformal invariance imposes severe 
constraints on f±{X) but does not fully determine them. To remedy this problem, two 
component amplitudes were computed by Feynman diagrams. We will see below that the 
Grassmannian integral is capable of reproducing the full 6-point superamplitude without 
any extra input. 

Cyclic gauge Let us begin by taking what we may call a "cyclic" gauge. 

^ C21 1 C23 C25 0^ 
C = C41 C43 1 C45 . (5.26) 
\C61 C63 C65 1/ 

The bosonic delta functions can be "integrated" by using the following relation, 

J d^^'^cS^Xr + CrsXs)S%cc^ + l)F(c) = (P) (F^) + F(c*_)) , (5.27) 

where are the two solutions of the delta-function constraints for the (3 x 3) matrix 
Crs' The barred (unbarred) indices label even (odd)-numbered particles, respectively. Both 
types of indices are defined modulo 6. Concretely, the solutions take the form 

. * X _ -(f|pi35|^) ±i{f + 2,f-2){s-2,s + 2) 

{c^hs - . (5.28) 

Here, we defined ^135 = pi -\- ps p^. Substituting these solutions to the fermionic delta 
function determines a±: 

SHC ■ 7j)l^^^ = S%Q)5^{a±) , ai = - ^ {epgfiPqWf ± i ^pApiWv) ■ (5-29) 



The minors in the denominator determine the bosonic functions /±(A), 

1 



± 



M1M2M3 



(^25 ^41 ^63 )± 



(5.30) 
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The component amplitudes can be easily read off from the superamplitude, for example, 



-1 

/yO* yO* yO* A ' ( yO* ^ 

V ^63 ^25 ^41 / + V ^63 ^25 ^41 > 



(<^(Ai)<^(A2)<^(A3)0(A4)<^(A5)</.(A6)) = . . . ' . + . . ^Tl. ^ ^ 

V^63^25^41>'+ l^63^25^41>'- 

(<^(Ai)</.(A2)^(A3)0(A4)^(A5)^(A6)) = + ."y^f- ^ ^^^^^^^ . (5.31) 

1^63^25 ^41 V<-63^25^41y'- 

In ref. [2], Aq^ and 746(/) were computed by Feynman diagrams. It was shown that the 
two component amplitudes and superconformal symmetry are sufficient to determine the 
full superamplitude. Even after judicious use of color-ordered Feynman rules and efforts 
to simplify the result, the final expressions for Aq^ and Ag^^ given in ref. [2] are somewhat 
lengthy. There, the cyclic symmetry of the amplitudes results from summing over diagrams 
related to each other by cyclic shifts. 

In contrast, our expressions for Aq^ and A^^f) take an extremely simple form. The 
component amplitudes contain two terms related by "A-parity", and each term exhibits 
manifest cyclic symmetry. Remarkably, our results exactly match those of ref. [2] evaluated 
at on-shell momenta. Technically, the comparison was done only numerically. But, since 
we are comparing two rational functions of the A-variables with quadratic constraints, 
numerical checks for sufficiently many A configurations amount to a complete proof. 

Factorization gauge Next, consider using the "factorization" gauge, 

/l ci4 ci5 cieA 



C 



(5.32) 



1 C24 C25 C26 
1 C34 C35 C36 J 

The computation is almost the same as in the cyclic gauge. The full superamplitude still 
takes the general form (5.25): 

The fermions a± are now given by (f = 1, 2, 3; 5 = 4, 5, 6) 

0^± = ~ 2(p)^3)2 i^PQriPQWr ± i epqr {pq) V^r) • (5-33) 

The matrices can be expressed as 

. * ^ {f\pi23\s)Ti{r + l,r + 2){s + l,s + 2) 

{C±)fs = 7- ' (^-2^) 

(Pl23)^ 

where cyclic identification is understood separately for the barred and unbarred indices. 
Because Mi = 1 in this gauge, the functions /±(A) take a simpler form 

[/±(A)]-i = MiMaMsl,^,^ = iTi)icU4,U 

_ .^ ((lbi23|4) T i(23)(56)) ((3|pi23|6) ^ i(12)(45)) 
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As a result, the component amplitudes also become simpler, for instance. 



A 



—I 



((2|pi23|6)-i(31)(45))3 



((l|pi23|4) - z(23)(56))((3|;)i23|6) - i(12)(45)) 
((2|pi23|6)+i(31)(45))3 1 



((1|P123|4) + z(23)(56))((3|pi23|6) + i(12)(45)) J ' 



(5.36) 



We have confirmed that the component amplitudes obtained in the two different gauges 
are indeed the same. The manifest appearance of the internal propagator 1/(^123)^ in 
eq.(5.36) suggests that the factorization gauge is more suitable for the study of the recursion 
relation. In fact, we will recover eq.(5.36) from the recursion relation in the next section. 
For a later convenience, we show two explicit examples in this gauge 



5.2.2 8-point amplitude 

Generality The counting (5.9) indicates that there is a one parameter family of solutions 
to delta function constraints. Following the general discussion above, if we fix the GL(4) 
gauge by filling up four of the columns of the C matrix by the identity matrix, the remaining 
(4 X 4) matrix Cfg defines (±z times) an element of 0(4, C). Consider the condition 



We can understand it geometrically as follows. We construct two sets of orthonormal basis 
for C^, one associated to (— z)Af (call it e^) and the other associated to (call it e^). 
We define ei and 62 such that ei is parallel to and 62 is parallel to A^ minus its 

projection onto ei. Then, 63 and 64 are defined to span the 2-plane orthogonal to ei^2- The 
other set {e^} is defined in a similar way with {—i)Xf replacing A^. 

The momentum conservation (5.38) requires that icfs map (61,62) onto (ei,e2). After 
imposing this condition, we need to map the 2-plane spanned by 63^4 to the 2-plane spanned 
by 63^4. There is precisely an 0(2) freedom in such a map, and that explains the origin of 
the single integration variable. 

The group 0(2) contains two copies of S0(2) related to each other by orientation 
reversal. Our experience with the integral formula in the 6-point example suggests that 
we should include both copies to find the correct answer. For a given fixed orientation, we 
parameterize the SO (2) rotation by r = e^^. Since we are working with a contour integral, 
9 and r should be understood as complex variables. 

Concretely, we can write the pair of one-parameter solutions c±{t) as follows. 



(<^(Ai)V'(A2)V;(A3)</'(A4)V5(A5)V'(A6)) 



(<^(Ai)(/.(A2)V;(A3)</'(A4)<^(A5)V(A6)) 




(5.37) 



(5.38) 



ic+(r) = eiej + 6263 + r(e+e^) + r ^(e-e^) , 
ic_(r) = eief + 6262 + r(e+e!^) + r~^(e_e!^) , 



(5.39) 
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where e± = ;^(e3zbze4), e± — -^{c'^^ie/^). The orthogonahty condition cc^ — —1 ensures 
that the minors Mi are at most quadratic in Cfs' Note that even for the minors quadratic 
in Cfs, due to the determinant structure, 

M{t) = Cf^si {r)cf^s2 (^) - Cfi52 (^)cf25i (t) , (5.40) 
the terms ah cancel out, so that the zeroes of the minors can be computed exphcitly: 



M(t) = aT-^ + /3 + 7T = =^ r= f^^VP^ 4^7 ^^^^^^ 

27 

The unpleasant square-root factors — 4^7 will all disappear once we sum over residues, 
and the final result will be a rational function of (ij) {i^j = 1, • • • 8). 

Since we reduce the matrix integral to an ordinary contour integral in the r-variable 
by solving delta function constraints, it is important to compute the Jacobian factors. As 
explained in appendix D, the Jacobian factor is fairly simple (po = J^rPf)'- 



V^27ri V/c+ JcJ r Mi{t)M2{t)M3{t)M^{t) ' 



27ri V/c+ Jc.J r Mi{t)M2{t)M3{t)M^{t 
where is given by 

= _^^^y \<^pqfs{pq)i4 + Crtvi){rii + Csurii)- (no sum in /) (5.43) 

The equahty in eq.(5.42) holds up to an overaU constant factor that cannot be determined 
by the Grassmannian integral. 

Factorization gauge Our earlier discussion of the 6-point amplitude shows that the 
amplitude in the factorization gauge is the simplest and most suitable for comparison with 
the recursion relation. For simplicity, we present the result only in the factorization gauge, 
(f=l, 2, 3, 4, 5 = 5,6,7,8) 



C 



/l Ci5 Ci6 Ci7 Ci8\ 

1 C25 C26 C27 C28 

1 C35 C36 C37 C38 

\0 1 C45 C46 C47 C48 / 



(5.44) 



The minors appearing in the denominator are given by Mi = 1 and 

M2;± = -ci5;±(r), M3;± = detP'J^J^^'^J^M , M4;± = ±C48;±(r), (5.45) 

\C25(r) C26(r) ' 



± 



with c±(r) defined in eq.(5.39). 
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For concreteness, we focus on a particular component amplitude A^^^^^^^^. This 
choice simplifies the evaluation of the contour integral considerably because the fermionic 
delta function in the numerator produces 



det 



C37 C38 
C47 C48 



3 

= (±M3±(t))^ , 



which removes poles corresponding to Ms{r). Note that the removal of the poles from 
Ms(r) is gauge independent. The gauge choice (5.44) further simplifies the calculation since 
Mi(t) = 1. We should stress that the net number of poles must be gauge independent; 
the poles for Mi(t) are merely pushed away to zero and infinity. 

The amplitude thus becomes 

A 1 1 // dT Mj^r) , / dr (-l)M|_(r) ^ 

A^^^,^,^^ - V M,^ir)M,^{r) + L V M,.{r)M,.ir) ) ' ^'"'^^ 

This looks much simpler than the generic integral (5.42). Let us introduce the notation 

{2}±, {4}±, {0}±, {(^}±, (5.47) 

where {0}±, {oc}± denote residues evaluated at r = and r = oc, {n}± is the sum of the 
two residues for the n-th minor. The sum of all residues should of course vanish, 

{2}+ + {4}+ + {0}+ + {^}+ = = {2}_ + {4}_ + {0}_ + {oo}_ . (5.48) 

To determine which linear combinations of the residues yield the correct 8-point am- 
plitude, we make use of discrete symmetries of the amplitude. First, consider a particular 
"A-parity" tt which acts as 

^ : Ai ^ -Ai . (5.49) 

Second, we consider the following permutation, 

a : {1,2,3,4,5,6,7,8} ^ {4,3,2,1,8,7,6,5}. (5.50) 

The amplitude at hand, A^^^^^^^^^ is invariant under both transformations; vr-invariance 
due to the fact that (^(Ai) is a boson, and cr-invariance due to the charge conjugation 
symmetry. On the other hand, the residues transform as 

vr : {2}± o {2}^, {4}± o {4}^, {0}± o {oo}^ , (5.51) 
a : {2}± o {4}±, {0}± o {0}±, {^}± o {^}± . (5.52) 

These transformation properties show that there is a unique combination of residues to 
give the correct amplitude, namely, 

Ah4>^^H^ = {0}+ + {0}- + {oo}+ + {oo}_ (5.53) 
= -({2}+ + {2}_ + {4}+ + {4}_). (5.54) 
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To emphasize the discrete symmetries, we can also write 



-^^m<l>H^ = (1 + vr) • ({0}+ + {0}_) = -(1 + vr) • (1 + a) • {2}+ . 



(5.55) 



After some algebra, we can write down the residues explicitly, 




(5.56) 



where n32±,?i34± are given by 



«32± = -V-^'i234((34)(6 b78|l) T (78)(2|p34|5)) 

-((12)(56) ± (34)(78))(l|p234|5) +p?234(6|P78i|2) , 

"34± = V-^'?234((56)(3|pi2|8) T (12)(4|p56|7)) 

-((12)(56) ± (34)(78))(4|p567|8) T J3?234(3|P456|7) . 



The numerator in eq.(5.57) has the expansion 



(^32+^34+ + ns2-n34-) = y -Pus^i' ' ' ) + (' ' ' )• 

The second term has eigenvalue ( — 1) under the vr-transformation, thus it is projected out 
by (1 + tt). To summarize, we obtain a remarkably simple final answer 



up to an overall constant. The overall factor 23^^678^234^567 eq.(5.57) suggests that 
the 8-point amplitude can be understood as a sum over off-shell 4- and 6-point sub-diagrams 
joined with an internal propagator, I/P1235 1/^678' 1/^234 1/^567- "^^^ absence of factors 
such as 1/^345 or 1/^455 also can be expected, since the corresponding decomposition to the 
sub-diagrams does not exist. In passing, we also note that under the deformation given in 
eq.(3.13), the large z behavior of ^^^^(f)^(f)^^{z) agrees not only with the Feynman diagram 
analysis in section 4, but also with the Grassmannian formula analysis in section 5.1. 

6 Examples of recursion relation 

In this section, we give a few examples of the recursion relation derived in section 3. We will 
match the recursion result with that derived in the previous sections. The checks are made 
up to an overall constant for the n-point superamplitude. In principle, by carefully defining 
the notion of color-ordered amplitudes, one could determine the overall coefficients without 
any ambiguity, both in the recursion relation and in the Feynman diagram analysis. As 
for the Grassmannian integral formula, there is no known way to determine the overall 
coefficients. We will not keep track of these coefficients in this paper. 



-(f)^^(f)i;(f)(f)i;(^l^^'2, A3, A4, A5, A6, A7, As) 



^^123^678^234^567 



K234 (-(2|P34|5)(3|P456|2>7|8) + {3\pi2\8){5\pe\p78l\2)) 



+ (12)(34) ((78)(l|p234|5)(4|p56|7) - (56)(8|2>i23|4)(6|j>78|l)) 

- (12)(56)2(l|p234|5)(8|pi2|3) + (34)(78)2(8|pi23|4)(2|;)34|5)}, (5.57) 
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Figure 7. Recursion relation applied to the 6-point amplitude. 



Example 1: (6) = (4)^(4) 

To be specific, consider a component amplitude (5.36) computed in the factorization gauge, 
which we reproduce here for the reader's convenience 

((2|pi23|6)-z(31)(45))3 



A 



(Pi23)2 L((l|pi23|4) -i(23)(56))((3|pi23|6) -i(12)(45)) 
((2|pi23|6)+i(31)(45))3 



((I|pi23|4) + i(23)(56))((3bi23|6) + i(12)(45)). " ^^'^^ 

To compare this with the recursion relation, we choose to deform (Ai, Ag). This deformation 
admits only one factorization channel (123:456) as demonstrated in Fig. 7. The recursion 
relation (3.26) reads 

A6(Al,A2,A3,A4,A5,A6) 

" (msF / ^^^(^^^l'^2)^4(Al(^D,A2,A3,A/)^4(iA/,A4,A5,A6(Zi*)) + (^l* ^4)) , 

(6.2) 

where H is the function defined in eq.(3.27) for even /. The component amplitude A^^^^^^ 
can be obtained from the superamplitude. Noting that fji{z)fjQ{z) — ryirye, we find 

^^^^^V^v^l-^i' A2, A3, A4, As, Ae) 
i 



bl23) 



(^H{zl,zl)A^^^^{\i{zl), A2, A3, \f)A^^^^{i\f, A4, As, A6(2:^)) + {zl 



(6.3) 



To obtain the explicit form of the amplitude, we need the expression for 4-point compo- 
nent amplitudes given in eq.(5.23). Due to the small number of external legs, the positions 
of poles take a simple form. 



(16)2 - ((23) - (45))2 



^2 



(16)2 - ((23) + (45))^ 



^ ((1| + i(6|)p45(|l) + i|6)) ' ((1| + i(6|)p45(|l) + m 

After repeated use of momentum conservation and Schouten's identity, we obtain 

H{^i^^2)A^^^^{zl)A^^^^{zl) = 



((2|pi23|6)+z(31)(45))2 



((l|pi23|4) +i(23)(56))((3|pi23|6) +i(12)(45)) 
((2bi23|6)-i(31)(45))3 
((lbi23|4) -i(23)(56))((3|pi23|6) -i(12)(45)) ' 



(6.4) 
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Figure 8. Recursion relation applied to the 8-point amplitude. 

Thus, we have verified that the amplitude obtained from the recursion relation (6.3) is the 
same as the amplitude computed directly from the Grassmannian integral (6.1). 

Obviously, we can check the recursion relation against the Grassmannian integral for 
many other component amplitudes and different factorization channels. For instance, if we 
deform (Ai, A3), there are two channels labeled by (561), (345). Schematically, the results 
can be expressed as 

^6^{^i, A2, A3, A4, A5, Ae) = A4^(A5, Ae, Ai, A/) A4^{iXf, A2, A3, A4) 

+ A4^(Xs, A4, A5, Xf)ik A4^{iXf, Ae, Ai, A2) , (6.5) 

^6(/)(Ai, A2, A3, A4, A5, Xq) = A4(p{X^, Ae, Ai, A/) ^ A4^{iXf, A2, A3, A4) 

+ ^4(/)(A3, A4, A5, A/) ^ A4^{iXf, Ae, Ai, A2) . (6.6) 

To check these relations in practice, care should be taken to assign the correct factors of 
(±z) on the right-hand side. 

Some component recursion relations can mix external particles. Let us show an explicit 
example. We are still deforming (Ai, A3). The recursion relation for A^^^^^^ can be read 
from the general recursion relation (3.19) by integrating out T]2viv5 of Aq. On the right 
hand side of the recursion relation, the component amplitude of A4 which depends on fji 
in turn depends on 773, thus it leads to 

^i,i;4>(P4>(p{^i^ ^2, A3, A4, A5, Ae) = -Ag\^3^3^3 

= s^{z*)A^^^^{X^, Ae, Ai, A/) ^ A^^^^{iXf, A2, A3, A4) 
+ c^(^*)^^./>V5V^('^5, Ae, Ai, Xf) ^ A^^^^{iXf, A2, A3, A4) + • • • . (6.7) 

where c(2:*), 5(2:*) are defined by c(z*) = — ^^i^*) = 2i — ' ^* zeros of 

the on-shell condition of the channel p'j^{z*) = 0. The first line of eq.(6.7) comes from the 
expansion of T]2Vivh while the second line from 772^3775. Note that in the first line, the 
external particles are changed from '^(Ai) • • • 4^{Xs) • • • to 0(Ai) • • • '^(^3) • • • . 

Example 2: (8) = (4) ^ (6) + (6) ^ (4) 

We apply the recursion relation to the 8-point component amplitude A^^^^^^^^ (5.57). If 
we choose to deform (Ai, Ag), only two factorization channels contribute to the recursion 
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relation as demonstrated in Fig. 8. Schematically, the recursion relation takes the form 

^4>^^P(P^P(P4>^p{^i^ A3, A4, A5, Ae, A7, As) 

= ^0V^V^(/)('^i' -^2, A3, A/) ^ A^^^^^^{i\f, A4, A5, Ae, A7, Ag) 
+ ^^v^v^(/)V^v^('^i5 -^2, A3, A4, A5, A/) ^ A^^^^(zA/, Ae, A7, Ag) . (6.8) 

To evaluate the right-hand side, we use the 4-point amplitudes in eq.(5.24) and the 6-point 
amplitudes in eq.(5.37). 

By numerically computing both sides of eq.(6.8) for a large number of {A^} configura- 
tions satisfying the overall momentum conservation, we have confirmed that the recursion 
relation holds. Since we are comparing two rational functions of complex variables, the 
numerical check amounts to a complete proof of the on-shell equivalence. 

It is worth noting that, if we express the 8-point amplitude as the sum of residues 
(5.54), the two channels of the recursion relation in eq.(6.8) separately match the residues 
in the following way, 

- ({4}+ + {4}-) = A^^^^(Ai, A2, A3, A/)^A^^^^^^(zA/, A4,A5, Ae, A7, As) , 

- ({2}+ + {2}-) = A^^^^^^(Ai,A2, A3, A4, A5, A/)^A^^^^(zA/, Ae, A7, As) . (6.9) 

One may expect it from the form of minors in the factorization gauge, M2 ^ C15, M4 ^ C4s, 
in eq.(5.45). 

7 Dual superconformal symmetry for all tree amplitudes 

In this section, we will use the recursion relation described in section 3 to demonstrate 
that the dual superconformal symmetry of the four-point amplitude can be extended to all 
on-shell tree- level amplitudes in the ABJM theory. As we will explain in section 7.2, this 
boils down to showing that the amplitudes are covariant under dual inversion. 

The proof is based on induction: assuming that An with n <2k transform covariantly 
under dual inversion, we demonstrate that each term in the recursive construction of A2k 
will transform in a universal covariant matter, and thus so will A2k- Since it was already 
shown that the four-point amplitude transforms covariantly under dual inversion [5], this 
proves covariance for all ABJM tree amplitudes. Similar proofs were given for the maximal 
super Yang-Mills theories in four-, six- and ten-dimensions [19-21]. 

In the next subsection, we will review the construction of the dual superspace for the 
ABJM theory, and explain how to perform inversion in this space. We present the proof 
of dual superconformal symmetry in section 7.2. 

7.1 Dual space and dual inversion properties 

In ref [2] , it was demonstrated that the four and six-point amplitudes of the ABJM theory 
have Yangian symmetry. Yangian symmetry of the scattering amplitudes was then shown 
to be equivalent to superconformal symmetry plus dual superconformal symmetry [5]. In 
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particular, the nontrivial level-one Yangian generators were matched with generators of 
dual superconformal symmetry, which are defined in a dual space parameterized by the co- 
ordinates (x, y, 0). The dual space coordinates are related to the on-shell space coordinates 
(A, 77) as follows: 



a/3 


a/3 


a/3 


a/3 
Pi - 




nia 




nia 










.ylJ - 

yi+1 - 


- J.IJ 
1 





(7.1) 

where X2k+i = ^i, = ^i? I/2/c+i = Vi- Note that the y coordinates are Grassmann- 

even. The dual space coordinates are defined such that supermomentum and R-symmetry 
are automatically conserved: 

5Zpi = 5Z9i = 5Zn = 0. (7.2) 



Inversion / acts on the dual space as 

a/3 a/3 nAanB/3 

The spinor indices (a, /3) are raised and lowered using the antisymmetric 2-index e tensor. 

As noted in ref. [5], the coordinates (xi^Oi) are sufficient to define (Xi^iji). Given the 
x's, one can determine A's using the first line in eq.(7.1) (up to sign ambiguities). These 
sign ambiguities also appear when translating from spinors to dual coordinates via 



(n + l) = ±^-x2.^2. (7.4) 
Given the 0's, one can then determine the 77's using the second line in eq.(7.1): 

r?f = (^^'^y^^^^+^P , (7.5) 

Finally, given the 77 coordinates, the y coordinates can be determined using the third line 
in eq.(7.1) up to some reference point yo, which the amplitudes do not depend on since 
they are invariant under translations in the dual superspace. Hence, the amplitudes can 
be parameterized using only x and 9 coordinates. Using this parametrization, only half of 
the dual supersymmetry is manifest. In principle, it should be possible to write amplitudes 
using all three coordinates (x, 0) in such a way that all of the dual supersymmetry is 
manifest, but for simplicity, we will only use x and coordinates to write the amplitudes. 

Note that the dual coordinates in eq.(7.1) are defined up to a constant shift. For later 
convenience, we fix this ambiguity by choosing^ 

X2 + X2k = 0, 6>i = qi. (7.6) 



^Note that this choice can always be made by taking the origin of the dual space to be at the midpoint 
of the line (x2,X2fc). 
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For this choice, 



Pi - P2k Pi + P2k Pi + P2k 

Xl , X2 -X2k , . . . , X2fc = , 

Oi = qi, ^2 = 0, . . . , 62k = ^1 + Q2k- 



(7.7) 



Furthermore, if we apply the deformation in eq.(3.13) to legs 1 and 2k, this only shifts the 
point (xi, 6^1) in the dual space: 



xi{z) , 9i{z) qi{z), 



Xi{z) = Xi, 9i{z) 



for z > 1. 



(7i 



Although Xl is deformed to its norm does not change, i.e. x^ [z) = xl- This IS one 

of the advantages of the choice in eq.(7.6). Another implication of this choice is 



x^ — X' 



2k 



(7.9) 



The inversion properties of (A, 77) are given by [5]: 



\/{Xi+lY{XiY ^J{Xi+lY{Xi)'■ 



{e^ = ±1) 



^1^1+1 



(7.10) 



where ej represents the sign ambiguity of the inversion rules, which is related to the sign 
ambiguity in eq.(7.4). The shifted dual coordinates obey these transformation rules if we 
choose X2 + X2k — and eie2k — —1: 



xfiz) 



where [x • 9)^ — x'^^9^. This can be demonstrated straightforwardly: 



I[xiiz)] 



2^2 



{xi ■ Xl - ieie2kX2k^X2kf_ -2 , (a^i • Ai + ieie2kX2k ■ Mk)"^ 2 



^X -^X ^2 

(Ai -iX2kf -2 , (Ai + iA2fc)^ 2 



^X-^ X'2 



4x1 



-z-^ + 



4x1 



xi{z) 



x^ 



(7.11) 



In obtaining this result, we noted that X2 = and used the explicit form of xi, X2, X2k in 
eq.(7.7). If we had chosen eie2k = 1, then eq.(7.11) would read 



I[xi(z)] 



xi{l/z) 



xt 



(7.12) 
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For simplicity, we will choose eie2k = —1- Choosing eie2k = 1 does not change our conclu- 
sions. 

Having defined our conventions for the dual space, we will now proceed to the proof 
of dual superconformal invariance. 

7.2 Proof of dual superconformal invariance 

We would like to show that when the on-shell amplitudes are written in terms of the dual 
superspace, they transform as follows under dual inversion: 

2k 
i=l 

or equivalently, 

2k 

n^?/2fc, A2k = f2kSHP)6%Q), (7.14) 

i=l 

where we use the identity I[6^{P)6^{Q)] = 6^{P)6^{Q) which is proven in ref. [5]. It is easy 
to verify eq.(7.14) for the case 2k = 4:^ 

(7-15) 

i=l 

Under a dual inversion-translation-inversion, eq.(7.13) implies that the superamplitude 
transforms as follows: 




(7.16) 



Hence, if we define the dual conformal boost generator to be 

^ 2k 

i=l 

then this will be a symmetry of the amplitudes. Furthermore, it can be shown that K 
matches a level- 1 Yangian generator and that all the other nontrivial dual superconformal 
generators can be obtained by commuting K with the ordinary superconformal generators 
(when acting on on-shell amplitudes) [5]. In summary, if we can prove that eq.(7.14) holds 
for all tree-level amplitudes, this implies that the amplitudes enjoy dual superconformal 
symmetry. 

^Note that there may be sign ambiguities when one combines square roots, since generically ^/A^/B = 
zLa/AB. This ambiguity can be removed by keeping track of the phase for each factor in the square root. 
In the end, terms with the same factors appearing in the square root wih have the same phase. 



HA2t] = 



\ 



i\hii = 




\ 
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Figure 9. The recursion in dual coordinates for 2k = 6. Note that the amphtude is deformed by 
shifting xi to xi{z). The shifts creates poles in the P123 channel and the residue is proportional to 
the product of four-point tree amplitudes at shifted kinematics. 



The proof is based on induction. Assuming that 



I[f2lixi)] = f2li^) 

xf 



21 



\ (7.18) 

\| i=l 

for all 2/, with I < k, we will use the recursion relation defined in section 3 to show that 
f2k inverts the same way. Recall the recursion relation given in eq.(3.26): 



where / labels the different channels and Al and Ar are lower-point amplitudes. We 
have Pf = Pj+i + . . . + P2k — ~{Pi + . . . + Pj) and {±2:^^, ±^2} are solutions of Pf{z)'^ = 
(Pi(^) + P2 + • • • + Pj)^ — 0- The function H(a,b) is given in eq.(3.27). Combining 
SHQl)S%Qr) = S^(Q2k)SHQR). one can extract /s^: 

d^r^\[5\QR)H{zl.j,zlj)fL{zlj,r^^ + ^ 4,/)]- (7.19) 

£ J Pf 



f 



We will now deduce the inversion property of /2/c by applying a dual inversion to 
eq.(7.19). As explained in section 7.1, the deformation corresponds to a shift in the dual 
coordinate xi. We illustrate this in Fig. 9 for the case 2k = 6. In terms of the shifted dual 
coordinate, 

/l(^*) = /l(xi(2:*),X2,...,xj+i), /i?(2:*) = /i^(xi(2:*),Xj+i,Xj+2,...,^2/c), (7.20) 

where we suppress the 0i coordinates for simplicity. By assumption, the functions /l and 
/r invert as 

/[/l(^i(^*), . . . , Xj+i)] = ^J {xKz'') . . . xj^^)fL{xi{z''), . . . , x^-+i), 

I[fR{xi{^n, . . . , X2k)] = ^(X2(z*)...x2,)/^(xi(z*), . . . , X2k)- (7.21) 
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The inversion of the propagator term in eq.(7.19) is given by 



A simple calculation shows that 



ge/jx(<9i - 6'j+i)^"A/^Q.((9i - 0j+iY^\f^p{9i - dj+i)^'^\f^^ . (7.23) 



To see how this expression inverts, we note that 



= ^=i==(^i - ei+i)'"Ay,« . (7.24) 

In the second hne, we used the identity Xi ■ Xi — xi+i ■ Aj with i — j + l,i + 1 — 1, Xi — \f. 
Hence, we find that 

/[ [ d^TjS^Qn)] = \ / d^TjS^'iQR). (7.25) 

xl{z*)x]^^^ xl{z*)x]^^ J 

The only remaining piece to invert in eq.(7.19) is H{z\, 2:2). From the observation that 

one can easily see that Pf{z)'^ = is equivalent to I[pf{z)'^] = 0. Thus 
so 

I[H{zlz;)] = Hizlz;). (7.27) 

Combining equations (7.19), (7.21), (7.22), (7.25) and (7.27), and using = x^, we 
conclude that 



2k\ 



2k 



A U^U2k. (7.28) 
\ i=i 



Hence, if inverts according to eq.(7.14) for all I < k, then so does f2k- Since /4 
satisfies this property, this completes the proof that all tree-level ABJM amplitudes are 
dual superconformal invariant. 
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Figure 10. A sample unitarity cut of four-point two- loop amplitude. Note that the dual points X2 
and xe, while adjacent, are not labeled successively. 



8 Dual superconformal symmetry of loop amplitudes 



In this section, we will demonstrate that the dual conformal properties of the tree-level am- 
plitudes can be extended to the cut-constructible parts of loop amplitudes via generalized 
unitarity methods [22-25]. To make a statement about the complete loop-level amplitudes, 
one needs to understand how the regulator modifies, or in some cases breaks, the symmetry. 
This is discussed in greater detail in the conclusion. 

We will follow the argument which was used to establish dual conformal symmetry for 
six-dimensional super Yang-Mills [20]. The proof involves showing that all non- vanishing 
unitarity cuts of a given loop amplitude, with cut propagators restored, invert in a universal 
fashion. Since all cuts invert in the same way, so must the amplitude, and this leads to the 
following loop-level statement: 



/[A 




i.i) 



where is understood to be the L-loop amplitude prior to integration, which is defined 
without a regulator. 

For loop amplitudes, the labeling of the dual space cannot always be done such that ad- 
jacent regions are labeled successively. For example, xq and X2 in Fig. 10 are non-successive 
yet they are adjacent regions. Thus we utilize the following more general notation: 



X, 



a(3 



a(3 _ a(3 

^3 ~^{ijy 



-i{ij} 



(8.2) 



In previous tree-level discussions, we had j — i^l. 

We begin by converting integrals over 77's to integrals over ^'s. The conversion is done 
using the following identity: 



JJ 5*^ {Or -9s- \{rs}'n{rs}) 



{rs} 
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where a runs over the tree diagrams in the cut, k runs over regions in the dual space, and 
{rs} runs over all cut lines and external lines. The number of regions for a loop amplitude 
is F = n + L, where L is the loop level. Since the translation from supermomentum to dual 
9 coordinates has an overall shift symmetry, the integration measure d^Ok is understood 
to include only F — 1 of them. To see the equality in eq.(8.3), one notes that the delta 
functions can be used to localize the 6(F — 1) integrals. Denoting the total number of 
lines, cut or external, as P, there are then 6(P — F + 1) delta functions left. For planar 
diagrams, P — V — F — 1^ where V is the number of vertices, or in our case, the number 
of tree amplitudes in the cut. Thus, we are left with a supermomentum delta function 
for each tree amplitude. In other words, by converting to the dual 9 representation, one 
automatically solves all supermomentum conservation constraints. 

Using eq.(8.3), one can rewrite the cut equation as: 



{ij} 



{ij} k 

xJ{S^er-es-\rs}ri{rs}) , (8-4) 
{rs} 

where {ij} runs over cut lines and {rs} runs over cut lines and external lines. After using 
the delta functions to eliminate the r^ijjj-dependence from each fa, the r/jjjj-dependence 
comes solely from the delta functions. The integral over then simplifies to 

J d\ij}S^ {Oi - Oj - \{ij}il{i3}) = {Oij\ij}) , (8.5) 

where dij^{ij} — {6% — 0j)"^a{ij}^ have suppressed the SU(3) R-index. Plugging 

this result into eq.(8.4) gives 



cut 



k a {ij} 

X n i^r -Os- X{rs}V{rs}) , (8.6) 
{rs} 

where {rs} now only runs over the external lines, i.e. there are n of them. Furthermore, 
one can pull out an overall supermomentum delta function, leaving (n— 1) delta functions, 
which completely saturate the 9 integrations over the external regions. We are finally left 
with 



cut 



\ k / {ij} a 
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where now k now runs over the loop regions, which are regions 5 and 6 in the example 
shown in Fig. 10. 

Let us consider the inversion weight of each term in eq.(8.7): 



4f- 



• For each loop region A:, the 9k measure contributes a factor (x^ 

• Each cut leg {ij} contributes which comes from 6^ 



• Each tree-level sub-amplitude contributes yYl^xf, where i runs over all regions ad- 
jacent to the tree. 

After restoring the cut propagators, which invert as 

J — T _ ^ J 



1 




1 




= I 











the resulting object obtains the following factor after a dual inversion: 



U-l 11(4 



.2 \3 



(8.9) 



where xi. are loop regions. 

As a concrete example, consider the diagram in Fig. 10. One has 



/ 


"^Fig.lO 








cut- 



f ^2 \ 3 / ^2^2^2^2 / ^2^2^2^2 / ^ 
{X^Xq) ^^ X'^X2XQX^^^ X'^X^X^X^^^ X 



2 ^2 ^2 ^2 



( o^2 ^2 ^2 ^2 ^2 ^2 ^2 ^2 \ 3 /2 
\^X2XqXqX^X^X^X^X^) ' 



cut 



Xa 



I rf2 



10 



cut 



Furthermore, there are four propagators in Fig. 10, which invert as follows: 



^.10) 



X 



15 



2^2 



1-^5 



X 



15 



X 



53 



2^2 



3-^5 



X 



53 



X 



56 



2^2 



X (\X 



6-^5 



X\ 



56 



^62 



2^2 



X Q 



2-^6 



^62 



Ul) 



Thus, when ^Fig.io 



cut 



is combined with the cut propagators, the resulting object has the 



following inversion weight: 



Yi2i rfi2 rfi2 rfi2 



12) 



which matches the result in eq.(8.9) 

There is one more ingredient missing, notably the loop integration measure J d^/^, 
which is written in terms of dual coordinates as J d^x^.. In three dimensions, the loop 
integration measure will provide an extra inversion weight of Ilil^^^J"^' which exactly 
cancels the extra weight of loop regions coming from the integrand. Hence, we find that all 
cuts, with the cut propagators and loop integration restored, invert the same way as the 
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tree amplitudes do. We conclude that the cut-constructible part of the L-loop amplitude 
inverts as 



\ \^=l / 

From the discussion in section 7, it follows that the cut-constructible parts of the loop-level 
amplitudes enjoy dual superconformal symmetry. 

9 Discussion 

In this paper, we constructed a BCFW-like recursion relation for three-dimensional Chern- 
Simons matter theories. This recursion relation involves shifting the momenta of two 
external particles. Unlike the usual BCFW approach, the shift is nonlinear in the complex 
deformation parameter z. Furthermore, the poles of the shifted amplitudes are computed 
by solving a quadratic equation, rather than a linear equation. 

Using background field methods and the Grassmannian integral formula, we have ar- 
gued that the superamplitudes of the AB JM theory vanish when the deformation parameter 
goes to infinity, which is required for the recursion relation to be applicable. We have explic- 
itly checked the recursion relation for six-point and eight-point tree-level amplitudes of the 
AB JM theory by reproducing the results of Feynman diagram calculations and the Grass- 
mannian formula. Finally, we have used this recursion relation to prove that all tree-level 
amplitudes of the ABJM theory enjoy dual superconformal symmetry. Using generalized 
unitarity methods [22-25], we have further extended the dual conformal symmetry to the 
cut constructible part of loop amplitudes. 

For = 4 sYM, the combination of ordinary and dual superconformal symmetries, 
along with invariance under helicity rescalings, uniquely fixes the amplitudes [34]. Given 
that all the tree-level amplitudes of the ABJM theory have dual superconformal symmetry, 
it is natural to ask whether the amplitudes of the ABJM theory can also be fixed by the 
dual and ordinary superconformal symmetries. If this turns out to be true, it would 
probably imply that the Grassmannian formula proposed in ref. [4] indeed generates all of 
the tree-level amplitudes. 

Another important implication of dual superconformal symmetry is that the ABJM 
theory may have an amplitude/ Wilson-loop duality. Since the one-loop correction to the 
scattering amplitudes vanishes, it would be useful to have an explicit calculation of the 
two-loop correction to the four-point scattering amplitude to compare with the Wilson 
loop computation [42]. 

Our discussion of dual superconformal symmetry at loop level is restricted to the cut- 
constructible parts of the amplitudes, or more precisely, prior to evaluating loop integrals. 
Due to infrared singularities, the integrals are ill-defined in the absence of a regulator. In- 
troducing a regulator will generally render the symmetries anomalous. On the other hand, 
if one can choose a regulator and modify the symmetry generator in such a way that the 
regulator becomes symmetry preserving, this implies that the symmetry is preserved at 
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the quantum level. This was done for = 4 sYM by considering the dual symmetry as 
five-dimensional, with the extra dimension giving rise to a massive regulator [43-45]. Al- 
ternatively, one can perform the unitarity cuts in a higher dimensional theory and perform 
dimensional reduction. This approach was applied to = 4 sYM, for example, by first 
demonstrating dual conformal symmetry in 6d maximal sYM and then taking the higher 
dimensional components of the loop momenta to be regulators for the four-dimensional 
integrals [20]. Unfortunately, Chern-Simons matter theories do not arise from dimensional 
reduction of any higher-dimensional theory, so it seems that the best approach is to in- 
troduce a mass deformation [1, 54]. A first step would be to see if the unitarity cuts of 
mass-deformed amplitudes preserve an "extended" dual conformal symmetry. 

It also would be interesting to investigate if the recursion relation we proposed is ap- 
plicable to three-dimensional theories other than the ABJM theory. The next candidate 
to consider is the BLG theory [46-48], which has maximal supersymmetry and a similar 
structure to the ABJM theory. Similarly, one could ask if this recursion relation is ap- 
plicable to theories with less supersymmetry. Although various component amplitudes in 
the ABJM theory have bad large-z behavior, the theory has enough supersymmetry to 
guarantee that all the superamplitudes have good large-z behavior. In other words, it is 
possible to relate all the badly-behaved component amplitudes of the ABJM theory to 
well-behaved component amplitudes using the supersymmetric Ward identities. This may 
not be possible for theories with less supersymmetry. It would therefore be interesting to 
determine the minimal amount of supersymmetry that is required to have good large-z 
behavior of the superamplitudes. 
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Appendix 

A Conventions 

We follow the conventions used in ref. [2]. The SL(2,R) metric is 

The spinor contraction is implemented as 

rXa^-V'aX", ep^A^^Ap, e"% = ^", e"%^ = 5^. (A.2) 

The vector notation is translated to the bi-spinor notation and vice versa through three 
dimensional gamma matrices, 

X^P = X^K)"^ = -\{cTnapX^^ (A.3) 



with 



„_ , -1 \ 1 / -1 \ 2 ^ 1 



We list some useful identities 

K)a/3Kr^ = -2r?'^% (A.5) 

(cr^)a/3(crA()75 = ea7e,9<5 + e/J^Ca^ , (A. 6) 

_|_ ^ac^bf^de _|_ ^bc^af^de _|_ ^ad^bf^ce _|_ ^bd^af^ce^ ^ (A. 7) 

^ _ ^ -e«/3^7 (A.8) 

^[al3] = - ^,9a = eo,^^^7' (A-9) 

where x'^ — x^^x^. 



x"^xp^ = -x^S^, (A. 10) 



B Feynman rules for ABJM 

We start with the Lagrangian for the ABJM theory [49-51]: 
C =C2 + Ccs + A + ^6 , 

Ccs =e^^^tr Q^^a^^A + ^5^^.A^^A - - "^gA^A^A^ , 

- 2igHr (v^^Va^b'/'"^) - ^5'^^ (V'a^'^'/'^^b) + 2igHr {^P <P^'4> b) , 
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Figure 11. Matter propagators 



p 



Figure 12. Gauge propagators 



where g — ^J^i^k^ (f)^ — i(ja — '0'^^ and 

Note that the first two terms in C4 break the R-symmetry from 50(8) to SU{4:) because 
they contain e^^^^. The bosonic potential term Cq is O (0^) and will not play an essential 
role in the analysis of this paper. 

The gauge group is U(A^) x U(A^), where Ajj, and are the associated gauge fields, 
and (f)^ and transform in the (N, N) representation. The fields can be expanded in 
terms of matrices as follows: 

where T", a = 1, ...,N'^ - 1, are hermitian generators of SU(A^) satisfying tr [t^T^^ — 5^^, 
T^" is l/\/]V times an TV x TV matrix, f « = -^T"" for a = 1, - 1, and f = 73^^'- 
Hence, at the level of matrix representations, the main difference between the gauge fields 
and the matter fields is that the matter fields are not hermitian. It should be emphasized 
that the matter fields carry two different indices, one for each U(A^). 

Plugging these expansions into the kinetic terms gives 

The propagators for the matter fields are depicted in Fig. 11. 

In the background field formalism, there are two kinds of gauge symmetry, notably 
gauge symmetry of the background and gauge symmetry of the fluctuations. We use the 
background gauge symmetry to impose the following conditions on the background gauge 
fields: 

g-A^ = g-i^ = 0. (B.l) 
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After gauge-fixing tiie background, we still have the gauge symmetry of the fluctuations 
(under which the background is inert). We can therefore use the Faddeev-Popov procedure 
to introduce the following gauge-fixing terms for the gauge-field fluctuations: 



tr 



SI ?2 



(B.2) 



where a and a represent the fluctuations of the gauge fields and (z = 1,2) are gauge- 
fixing parameters. Note that these gauge-fixing terms don't preserve the background gauge 
symmetry. On the other hand, the background gauge symmetry is already broken by eq. 
(B.l), so there's no need to choose gauge-fixing terms which preserve the background gauge 
symmetry. 

The propagators for the gauge field fluctuations are then given by 

+ (B.3) 

where dz refers to the a/ a fields and a, h in the above equation are adjoint indices. Taking 
ii — Q gives Landau gauge propagators: 

The resulting propagators are depicted in Fig. 12. 

The color-ordered Feynman rules for the cubic and quartic vertices are depicted in Figs. 
13 and 14. They can easily be adapted to the background field approach by choosing at least 
two legs in each diagram to be fluctuation fields. The symmetry factors are unchanged. 
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Figure 14. Quartic vertices 



C Alternative forms of recursion relation 
Linear multi-line shift 

Consider four external particles labeled i, j, fc, I and define 

Ma = (A^ + iXj)^ , Pa = (Afc + iA/)^ . 

Also define fl = Xi — iXj and p = Xk — iXi. Note that P(^aUi3) — {Pi ~^Pj)ai3 P(aPp) — 
{Pk +Pl)a^- Now define the shift 

Xi^ Xi + -zp, Xj Xj + -zp, Xk^ Xk - -zp, Xi^ Xi + -zp. 
Under this deformation, the momenta shift as follows 

Pi + Pj -^Pi+ Pj + zq, Pk+Pl ^ Pk+Pl- zq, 
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where qa/3 = l^{aPp)' Note that pi pj -\- Pk -\- Pi is invariant under the shift. Also note 
that each momentum remains on-shell after the shift. It is straight forward to define an 
analogous shift of the fermionic coordinates. 

Near each pole of the amplitude, the amplitude factorizes into a product of two on- 
shell amplitudes connected by a propagator. Suppose that ij appear on one side of the 
propagator and k,l appear on the other side. In this case, the momentum in the propagator 
shifts linearly as follows: 

Pf ^Pf{z) =pf + zq. 

The poles of this propagator are therefore obtained by solving a quadratic equation. Now 
suppose that only particle i appears on one side of the propagator. Then 

Pf pf{z) =Pf + ZQa + z^qt , 

where (qa)a/3 = and {qt)^^ = IPaP^- Since q^ = q^ • q^ = 0, Pf{z)'^ is once again 

quadratic in z. 

Following the argument in section 3.1, it is not difficult to show that the recursion 
formula for the four-line shift is 

A{z = 0) = J2 [ ^'^^ (^(4' zJ)AL{zpAR{zp + (4 ^ zj 

where X]/ ^^^^ possible channels, z^ are the zeroes of p'f{z), and 

H{x,y) ^ 



x-y 



Note the similarity to the recursion formula for the two-line shift proposed in section 3.1. 
Although the recursion formula for the four-line shift looks a bit simpler, the two-line 
shift still has several advantages. For example, the recursion formula for the two-line shift 
involves fewer channels than the recursion formula for the four-line shift. 



Superconformally covariant form 

The BCFW recursion relation for = 4 sYM is known to take the simplest form in the 
super-twistor notation [52, 53]. By doing a similar manipulation in our case, we will show 
(rather heuristically) that the recursion relation derived in section 3 can be written in a 
superconformally covariant form. 

Review of the 4d result Let us review the ideas of refs. [52, 53]. They begin with 
the usual BCFW recursion relation in terms of amplitudes with the momentum conserving 
delta functions removed. Schematically, 
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Here is the value of the BCFW deformation parameter z satisfying [P/^(z*)]^ = 0. In 
translating eq.(C.l) into the twistor language, the first step is to insert the momentum 
conserving delta functions on both sides using the identity 

S\P) = 5HPl{z*) + Pr{z*)) = j dSS\PL{z*)-v)5\PR{z*) +p) . (C.2) 

Using the notation A — A5^{P)^ we have 

^-T.j'^'P {^L\pi{z*); -p]j2^R\Pj{z*y, . (C.3) 

The next step could be called "undoing the BCFW residue calculus" as it recovers the 
z-integral. Technically, we use the identity 

^ = -sgn([i|P,.|i)) I ^S{z[i\PL\j) + Pi) ^ -sgn{[i\PL\j)) J '^5{PL{zf),{CA) 
and use the delta functions to rewrite the formula such that 

■^=-E/t/ A^(p')sgn([ib|j))^L[p.(z);-p]^K[p,(z);+p]. (C.5) 

Reinstating the z-integral highlights the on-shell nature of the BCFW recursion relation 
in two ways. First, it removes the off-shell propagator (1/P^). Second, it restricts the 
integration domain of p to on-shell via 5{p^). The on-shell momentum integral J d^p5{p^) 
transforms into the J d?\d?\ integral at a later stage, which in turn translates into the 
twistor language. 

In eq.(C.5), the z-integral looks exactly the same as the BCFW contour integral, 
although with a different integration contour. It is not a coincidence; it follows from the 
identity, 

which holds at least formally for an arbitrary function g{z). 

Generalization to 3d It should be clear from the above review that the first two steps 
— inserting momentum conserving delta functions and undoing the residue calculus — can 
be straightforwardly carried over to three dimensions. We find 

= E / ^ / dhKp'')Xij{p)AMz); -p\Ar\p,{z)-+p] , (C.7) 
^ J z J 

with some function Xij{p) analogous to sgn([i|p|j)) in eq.(C.4). It must be a dimensionless 
function of p, Aj and Aj, but its precise form is not of interest in this paper. 

The next step is to make the replacement 

j (fp5{p^) j d^\, (C.8) 
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which follows easily from a change of variable such as 

pii = Ai Ai , P22 = A2A2 , P12 = 5 Ai A2 . (C.9) 

Supplementing the A-integral with the 77-integral, we can write down the superconformally 
covariant version of the three-dimensional recursion relation: 

= E/i^/ d^^''AXij{A)AL[Ai{zy,A]AR[Aj{zy,iA] 

= Y, [ d^^^AX,,{A) I -^z^^^ {AL[Af,A]AR[Af,iA]) . (C.IO) 

On the second line, Mij denotes the SO (2) deformation generator written in terms of the 
A- variables: 



M« = i(A«_A»). (C.U) 



In four dimensions, after similar rearrangements and half Fourier transformations, all 
the z dependence are absorbed into 

where Wi and Zj are four dimensional twistor variables for the legs i and j, respectively. 
It is not clear how to further simplify the z-dependence in eq.(C.lO). 

D Computational details 
Large-z behavior 

We demonstrate the cancelation of bad large z behavior for the four-point diagrams in Fig. 2. 
Note that the hard scalar and spinor have the same SU(4) R-index in these diagram. As 
a result, the 4-point contact term involving ^^^^^ does not contribute. 

The diagram in Fig. 2 (a) is given by 



In the large-z limit, pi z^q, tp{pi) zXq, and the above formula simplifies to 



q • k2 

Using eq.(A.7), it is not difficult to show that 

{q\ cj^^{ki){e^,pq^e^)^{k2) = {q\ ^{ki){k2 • qMk2) . (D.l) 



A similar gluing of amplitudes without the deformation was considered in ref. [2]. Strictly speaking, 
the superconformal covariance of this A-integral representation relies on the precise form of Xij (A) . 
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Hence, we get 

zg^{q\^{kiMk2) + 0{l/z) . 
The diagram in Fig. 2(b) is simply given by 

In the large-z hmit, this becomes 

Thus we see that in the large-z hmit, the terms of 0{z) cancel out. 
The 8-point calculation 

We present some of the complicated details in calculating the 8-point amplitude. We begin 
with a closer look at the orthonormal basis {e, e} defined in section 5.2.2. The possible 
ambiguities in the orientation of the basis can be lifted by choosing the signs of Si,£2,£3 
in the following identities 

Its') — Its') 

y-v^ V-Po 

^ V-Po ^ V-Po 

f \ f \ ^ f \ f \ e , {r\P0\s) f- ^ f- ^ , /- N /- N r {r\Po\s) 

[es)r[es)s + [e4)r[e4)s = Ors H 2 ' \^Vr[es)s + [e4)r[e4)s = Ors 2 ' 

^0 ^0 

(ei)r(ei)s + {e2)r{e2)s = ^ ^^^^2^"^^ (PO = ^Pr)' (D.2) 



Po 



r 



where po = P1234 for the factorization gauge, and po = P2468 for the cyclic gauge. In this 
paper, we always work with the choice £1 = £2 — £3 — 

The constraint equations imposed by the bosonic delta functions are 

One particular solution for the equations is 

c = -z(eief + 626^ + 6363 + 6464 ). (D.4) 

One can check that eq.(D.4) indeed solves eq.(D.3) using the completeness and orthonor- 
mality of the basis and the momentum conservation condition of A's. 

We use the following pair of one-parameter family of solutions {c+(t), c_(t)} (t = e^^) 

= «feef + ft4 + «€«') + e-AW'-). ( '^i'i) = ( 

\ e4[0) I \ — smt^ cosd' / \ 64 / 

\ 64(6^) / \ — sint' cost' / \ —64 / 
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To reduce the Grassmannian integral (5.42) down to an ordinary contour integral in r, we 
introduce the following new variables (4x4 matrix) {(jj)fi. 

{c)fp = • e)rp, 4x4 matrices {e)ip = (e^)^, {cj)ri = {cjf)i- (D-5) 

In terms of these variables, the solutions are given by 

(jj±{0) = {ei, e2, cos 0es + sinOe^, ^ sinOes ± cos Oe^}. (D.6) 

Ignoring the fermionic delta function for a while, the integral gets transformed to 

^'^^^ ^ J " M1M2M3M4 ' 

= ^y^" M,M,M,M, ' ^ = det(ei,e2,e3,e4) = ±l. 

= / rf^(cj)f3^^(^)f4 7.^ 7.^ 7.^ 7.^ ^ integrating out 8 variables (cj)fi, (cj)f2. 

(Po) J M1M2M3M4 

Now let the four barred indices be 

where {fi, f2, f3, f4} = {1,2,3,4} for the factorization gauge. We now integrate out six 
variables ^^^23, ^f245 ^fsS? ^f345 ^f43 5 ^f44- Then the bosonic integral becomes 

^'^^^ ~ ^K-i)3rfK-.)4 M,M2M3M4nK=2,3,4 2(-.,3-..4-C..,4C...3) ' ^""''^ 

Using the identity 

5(a;^=2 • t<^f3)^(^f3 • ^f^K^f^ • ^f-i) = <^^(-P)(^2r3)(r3f4)(f4f2), (D.9) 
the integral is further simphfied to 



(Po)^ J ^'"'■''^^'"'■'^^MiM2M3M4 ni^=2,3,4 2K-i3^^fK4 - C^fi4C^f^3) 
{pIY J dO 2a;f:i4 Ml M2M3M4 17^=2,3,4 2(^^ri3^^fK4 -^^ri4^^fK3) 



J Ml 

/ M1M2M3M4 ^ J ^^MiM2A 



1 I 

U}=Ul'l 



M2M3M4 ni^=2,3,4 2(Wfi3<^fK4 ~ Wfj4C<;f^3 

1 



M1M2M3M4 nK=2,3,4 2(cUf=i3Wf=^4 - CUf=i4W?^3) 

leyz^ [J Ml M2M3M4 ' J M1M2M3M4' 



(D.IO) 



In the first three lines, we used the shorthand notation (rr)^ = {r2rs){rsr4){r4r2). In 
the last line, we used a;fi3^fK4 - (^ri4(^rK^\uj=ojl = ±((^3)^ (e4)fK " (^4)fi (63)7^^) and the 
identity (D.2). Including the fermionic delta function and changing the integral variable by 
T = e^^, we obtain the final expression in eq.(5.42) by choosing some appropriate contours. 
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